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, ^ , A simple analytical expression for the electric dipole polarizability of the three- hadron 

bound system having only one stable bound state has been derived neglecting by 
the higher orbital components of the off-shell three-body transition matrix at the 
' energy of the bound state. As a case in point, we have estimated the electric dipole 

CO \ polarizability of the triton, using a cluster triton wave function and the Hulthen 

potential to describe the related p — n and n — d bound states. 
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1. Introduction 



1> \ Data on the electric polarizabilities of the lightest nuclei oe contain a valuable 



information on the nuclear force between nucleons. 
■ For the deuteron, the currently available values of q;£;(^H) obtained by the direct 

measurement of deviation from Rutherford scattering of the deuteron on a heavy 
2 ! nucleus below the Coulomb barrier^ and extracted from photoabsorption data^, 



a^(2H) = 0.70 ± 0.05 fm^ (Ref. 1) and ^^^(^H) = 0.61 ± 0.04 fm^ (Ref. 2) , (1) 

^ ■ are a little distinguished between themselves, 

'jjj . For the nucleus ^He, contrastingly, the corresponding values of the polarizability, 

5^ I obtained by the direct way^ and from experimental photoabsorption data^, 

aECHe) = 0.25 ± 0.04 fm^ (Ref. 3) and ^^(^He) = 0.130 ± 0.013 fm^ (Ref. 4) , 

(2) 

are distinguished by a factor of two. 

For the nucleus no measurement of the electric dipole polarizability has been 
performed to the present time. 

The polarizability of the nucleus "^He was experimentally found^"^ to be less by 
about an order of magnitude than that of the deuteron. 

Calculations of the deuteron electric polarizability, carried out with the re- 
alistic nucleon-nucleon interaction potentials, lead to the values of a£;(^H) = 
0.6328(17) fm^ (Ref. 7), being closer to the data of Ref. 2 in (1). 

Furthermore, examining theoretically the anisotropy of the deuteron deformation 
in the electric field caused by the tensor character of the n — p interaction, the 
separate components of the deuteron electric polarizability a'^^^', the longitudinal 
component (with the deuteron spin along the electric field) and the transverse 
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one a%, have been calculated in Ref. 8. (The above electric polar izability q;e(^H) 
is the averaged value of the components a'^', aE(^H) = |q;e + .) 

Computations of the scalar and tensor deuteron polarizabilities have also been 
performed in the framework of the effective field theory that uses space-time and 
global chiral symmetries of the quantum chromodynamics consistently describing 
pion propagation and relativistic effects^"^^. The results for the electric deuteron 
polarizabilities obtained in the cited works in the leading and next-to-leading 
orders agree with the values calculated in the traditional nuclear physics with the 
application of the potential models. 

The results of calculations of the '^He polarizability, a£;("'He) = 0.145 fm^ 
(Ref. 12), a^^'^Hc) = 0.153(15) fm'^ (Ref. 6) and a£(-^He) = 0.149(5) fm^ 
(Ref. 13), support the experimental result in (2) that has been obtained from 
the photoabsorption data (Ref. 4). Also, the results of calculations of the triton 
polarizability, q;e(^H) = 0.139(2) fm^ (Refs. 13 and 14), turned out to be close to 
those of the ^He polarizability^'^^-^'^. (A little larger value of Q;£;(^He) may be assigned 
to the repulsive Coulomb interaction between the two protons in the nucleus ^He 
causing the charge symmetry violation.) 

In the previous paper^^, we have predicted the value of the electric dipole 
polarizability of the only three-body lambda hypernucleus — the lambda hypertriton 
^H. It was found that aE^Ji) is close to 3 fm^ exceeding the polarizability of the 
ordinary three-body nuclei by an order of magnitude and even the recently measured 
polarizability of the nucleus ^He, ^^(^He) = 1.99(40) fm^ (Ref. 6). 

This paper is devoted to development of a method of direct determination of the 
electric dipole polarizability of the three-particle bound system, leaning upon solu- 
tion of the three-body problem at the negative bound-state energy without necessity 
of finding the continuum wave functions. In Sec. 2, a general formalism for deter- 
mining the polarizability of the three-hadron complex is worked out. The electric 
polarizability of the three-hadron nucleus is expressed in terms of partial derivatives 
of the bound-state wave function in momentum space and selected higher partial 
components of the three-body off-shell transition matrix. In Sec. 3, the developed 
formalism is applied in the case of a simple physically justified (cluster) model of 
the wave function of the triton (a bound complex of one proton and two neutrons). 
Neglecting higher partial components of the transition matrix we obtain a simple 
formula for estimation of the triton electric dipole polarizabihty. The results of 
corresponding calculations of q;£:(^H) based on the known low-energy data for the 
p — n and d — n interactions are described in Sec. 4. Conclusions are drawn in Sec. 5. 

2. General formalism 

Previously, on the basis of the three-body formalism of the effective interaction of a 
charged particle and a bound complex^^"^^, we have derived an expression for the 
polarization potential of the two-hadron bound complex that consists of charged and 
neutral hadrons (for example, the deuteron) starting immediately from the 

Faddeev integral equations^^. In the case that the interaction between the proton 
and the neutron composing the deuteron is central and 5'-wave, the electric dipole 
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polarizability of the deuteron is given by 




(3) 



where Cp is the charge of the proton, jipn is the proton-neutron reduced mass, 
IXpn — mpmn/rripn, rupn — mp + m„ {rup and m„ are the proton and neutron 
masses), V'd(^) = d%l)ci{k) / dk is the first derivative of the deuteron wave function 
in the momentum space in the variable of the relative momentum A;, and is the 
binding energy of the deuteron. The formula (3) is in agreement with the expressions 
for aE obtained in the case of the separable S-wave pair potential in Refs. 17, 20 
and 21. 

According to the expression for the electric polarizability of the two-particle 
bound complex (3), obtained assuming that the interaction in the P-wave orbital 
state is absent, the quantity a^; is completely determined by the wave function 
of the bound state of the complex. Hence, for the different interaction potentials 
producing identical bound-state wave functions (as an example, in the case of the 
two-body problem with the S'-wave local Hulthen interaction potential and with 
the S'-wave separable potential having the Yukawa formfactor), even if distinct the 
corresponding scattering wave functions, the expressions for the polarizability aE 
should be the same. 

In particular, from this it follows that a closer determination of q;e is anticipated 
in the event of fitting the potential parameters to bound-state data rather than 
scattering-state ones corresponding to higher energies. 

Here we obtain the expression for the electric dipole polarizability of the three- 
hadron bound complex in the case when the system can form, apart from the 
continuum, only one bound state. We start from the general expression for the 
polarizability of the complex , considering the low-energy scattering of the three- 
body bound complex (with the binding energy Bq) by the field of a particle 1 having 
the electric charge Ci. For simplicity sake assume that the complex consists of one 
charged particle 2 and two neutral particles (3 and 4). The initial kinetic energy of 
the relative motion of the particle 1 and complex is taken to be far lower than the 
breakup threshold energy of the complex. The effective potential of the interaction 
between the charged particle 1 and complex was found within the framework of the 
rigorous Watson- Feshbach formalism^^'^^ (for more details concerning application of 
the above formalism for the three-body system see Refs. 17 - 19, 21 ). The electric 
dipole polarizability aE was determined as the strength of the polarization potential 
at asymptotically large distances between the particle 1 and complex, pi, greatly 
exceeding the size of the complex. 




(4) 



where is given as 



= -2 < ^0 I (D2 • Pi)G'^(-Po)(D2 • Pi) I *o > , 



(5) 



Here \E'o) D2 = e2r2, and G'^{—Bq) are the wave function of the ground bound state 
of the three-hadron complex corresponding to the binding energy Bq (normalized 
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to unit, < ^0 I ^0 >= l)i the operator of the dipole moment of the charged 
particle 2 having the charge 62, and the truncated Green's operator of the complex 
G^{E) = QG{E) at the energy E — —Bq, respectively. The full Green's operator 
of the complex G{E) is given by G{E) — {E — Hq — V)~^, where Hq is the kinetic 
energy operator and V is the total interaction potential, V = V23 + f 24 + ^"34, Vij is 
the potential of the pair interaction between the particles i and j (the potentials 
are assumed to be energy- independent), Q = 1 — P, P is the projection operator 
onto the complex ground state, P =| >< ^0 |- The quantity Pi in Eqs. (4) and 
(5) is the radius vector specifying the relative position of the centre of mass of the 
complex (composed of the particles 2, 3 and 4) with respect to the charged particle 
1, = Pi/pi is the unit vector along p^, giving the direction of the external electric 
field, created by the particle 1, r2 is the radius vector of the charged constituent of 
the complex, the particle 2, relative to the centre of mass of the complex. 

Expressing the full Green's operator G{E) — [E — Hq — V)~^ through the 
transition matrix T{E), 

G{E) = Go{E) + G,{E)T{E)Go{E), (6) 

we write the operator G^{E), which appears in the expression (5) for the electric 
dipole polarizability of the three-body bound complex in the form 

(?«(£;) = (1 - P)[Gq{E) + Gq{E)T{E)Go{E)] , (7) 

where Gq{E) = {E — Ho)~^ is the free Green's operator. The transition matrix T{E) 
is defined by the Lippmann-Schwinger integral equation 

T{E) = V + VGo{E)T{E) . (8) 

In the case under consideration that the complex has only one bound state with 
the energy E = —Bo, the transition matrix may be written as the sum of the pole, 
T{E), and smooth, T{E), parts, 

T{E)^f{E)+f{E) , (9) 

where 

f{E) = , (10) 

I r„) = G„-'(-B„) I So) = V I *„> . (11) 
At the point E — —Bq the operator T{E) may be shown to have the form 

f(-i?o) = -ci |ro)(ro|+--- , (12) 

with 

ci = (^0 I Goi-Bo) I ^0) . (13) 

The factored term in (12) is the smooth part of the dominant partial component 
of the three-body transition matrix T{E) with zero relative orbital momenta. 
The additional terms in (12) indicated by the eUipsis contain the higher orbital 
components of T{E) ai E — —Bq, which are wholly smooth functions of E, however, 
being less important, they are disregarded here. 
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In view of (7) and (9), the operator G^{E) may be written as 

G^{E) = G^{E) + G^{E) + G^{E) (14) 
while the corresponding terms are given by 

G'^{E) = {l-P)Go{E), 

GQ{E) = {l-P)Go{E)f{E)Go{E), (15) 
GQ{E) = {l-P)Go{E)f{E)Go{E). 

The value of the matrix element 

(*o I (D2 • Pi)G«(£;)(D2 • pi) I ^'o) = (*o I (D2 ■ p,)G^{E){B2 ■ p,) \ ^o) 

+ (vl/o|(D2-Pi)G'«(E)(D2-Pi)|^o) 
+ (*o| (D2-Pi)G'Q(£;)(D2-pi) l^o) 

(16) 

appearing in the expression (5) at the energy E — —Bq is obtained with the use of 
the relations (10) — (12) in the second and third terms on the right-hand side of Eq. 
(16) and with evaluating the indeterminate form of the type 0/0 in the second term. 
Furthermore, it is easy to verify that the third term vanishes at the point E = —Bq. 
As a result, the expression for the electric dipole polarizability (5) becomes then 

aE = - 2{(^o I (D2-Pi)G'o(-Bo)(D2-Pi) I ^o) 

+ (^'o I (D2 ■ Pi) I ^0) ■ [ci(^o I (D2 ■ Pi) I ^0) 

- I (D2 ■ Pi)Go(-i?o) I ^0) - (^0 I Go(-i?o)(D2 ■ Pi) I ^0)]}. 

(17) 

Further simplification of the expression (17) for occurs if the wave function 
of the bound complex ^0 is characterized by a definite parity (as a consequence of 
the spatial reflection invariance of the interaction potential V). In such a case the 
matrix elements (^0 I (D2 ■ Pi) | ^0) , (^0 I Go{-Bo){'D2 ■ Pi) | and 
(^0 I (D2 • Pi)Go{—Bo) I ^0) in (17) vanish due to integration over the angular 
variables. The expression (17) then reduces to the simple form 

aE = -2(*o I (D2 • Pi)G'o(-So)(D2 • pj | *o) (18) 

containing only the free Green's operator at the negative energy E — —Bq. 

The bound complex being considered below is a three-hadron nucleus composed 
of the proton p (the charged particle 2), the neutron n (the neutral particle 3) and 
the neutral hadron h (the particle 4). The symbol h stands for n (the neutron) in 
the case of the triton and for A (the lambda hyperon) in the case of the lambda 
hypertriton. 

In the momentum representation the operators of the dipole moment D2 and the 
free propagator Go{—Bq) in the Eq. (18) are given by 



{kp\Goi-Bo)\k'p') = -i2.mk-k')5{p-p'){Bo + ^^+^y . 

where k and p are the Jacobi momentum variables describing the relative motion 
of particles p and n and that of the particle h with respect to the centre of mass of 

TT^nkp nipkn ^ 7T7./i(kp + k„) TTip^k/j (20) 
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Here kj is the momentum of the particle i, fipn,h is the reduced mass of the system 
(p, n) with the mass rripn = nip + m„ and the hyperon h with the mass mh, 
l-i'pn,h — TTipnmh/mpnh, nipnh = mp + TTin + mh, the binding energy of the three-hadron 
bound complex Bq is equal to the sum of the dcuteron binding energy Bd = K^/^Hpn 
and the separation energy of the hyperon h, Bh = K'fj2fid,hi = Bd + Bh, 
IJ'd,h = nidmh/mdh, md is the deuteron mass, rudh = rud + mh, \I'o(k, p) is the 

d_ 
9k 



normalized wave function of the three-hadron nucleus in momentum space, Vk ~ ^ 
is the gradient operator. 

When the full wave function in the general formula (18) has the form of the 

product of spatial and spin functions, we may omit spin variables from consideration. 
In the momentum space, substituting the expressions (19) into (18), the formula for 
the electric dipole polarizability of the three-hadron bound system may be written 
as 

r dkHn I (kp I P, • f^Vk + ^^Vp) ^pnh) P 
""^^ ^ h'c^ J (27r)6 ^ + ^ + So ^ ^ 

where "ifpnh is the spatial wave function of the three-hadron nucleus. 

We are reminded that the formula for the electric polarizability of the three- 
hadron bound system (21) has been derived assuming that the higher orbital 
components of the off-shell three-body transition matrix at the negative energy 
E = —Bq are negligibly small. Similar to the expression (3) for the electric dipole 
polarizability of the two-body complex, the obtained expression for the polarizability 
of the three-body complex is essentially determined by the first derivatives of its 
bound-state wave function with respect to the Jacobi momentum variables. 

Below, using a simple model wave function we employ our formula (21) to 
evaluate the electric dipole polarizability of the simplest three-body nucleus, the 
triton containing only one charged particle (the proton). 

It is worth noting that the parametrized analytical form for the Faddeev triton 
wave function generated with the Reid soft-core potential in Ref. 25 is quite suitable 
to use in the expression of the type (21) with the aim to calculate the triton electric 
polarizability of the triton. 



3. Model wave function of the triton 

In this article we study the electric dipole polarizability of the triton nucleus 
containing the proton (the particle 2) and two neutrons (the particles 3 and 4). 
To estimate the polarizability, we use the clustered (deuteron -|- neutron) triton 
wave function that must be antisymmetrized with respect to the identical fermions. 

For simplicity sake, without introducing the isobaric formalism, we have to 
antisymmetrize only with respect to the two neutrons. This antisymmetrization 
naturally occurs, if we shall restrict our consideration only the dominant symmetric 
part of the space wave function of the triton, since its spin wave function Xt = Xsms 
(with 5" = 1/2) is already antisymmetric relative to the permutation of the neutrons, 

XI 1 (3, 4; 2) = (a3/34 - a^Ps) «2. (22) 

where ctj and /3i are the spin functions of the particle i with the spin projections 
+1/2 and -1/2. 
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Thus, in the framework of the cluster model of the triton as a n + li system, the 
space part of the triton wave function must be symmetric in the interchange of the 
two neutrons, 



(23) 



where is the normalized bound-state wave function describing relative motion 

of the neutron and centre of mass of the deuteron, the Jacobi momentum variables 
= k and p4 = p are defined by the relations (20), 



m 



£■24 



p 1 1 ''pnn 



p, P3 = k - 



-p , 



(24) 



and the coefficient ensures normalization of the whole triton wave function to 1, 

7V„ = (1 + {Mk)Mp) + Mk24)MP3))y'^' ■ 

Substituting the model wave function (23) into the formula (21), we get the 
following expressions for the electric dipole polarizabilities of the triton. 



with 



q;b(^H) — a%{pnn) + a'^£{pnn) 



3 h m. 



pn 



p^dp 
2^ 



rrir 



Inip) 



pnn y 



d(l)n{p) 

dp 



(25) 



and 



+ 



Jn{p) [<Pn{p)Y 



pn J 



dkdp 



h m. 



pn 



r f dka 
J J 72^ 



(27r)6 + [C, 



m 



Mk)Mk24) 1 (Pl • P) [Pl ■ P3) 



pnn ^ 



+ 



\ 2 



'pn , 



dk dk- 



dp dp3 



0n(p)0n(P3)(Pl • ^)(Pl • ^24) 



24 



(26) 



(27) 



Here the hat signifies a unit vector, the functions Inip) and Jn(p) are defined by 

dke [Mk)? J ,„^ _ r dkk"^ [dMk)/dkf 



In(p) 



27r2 fc2 + [C„(p)]2 

[Cn{p)r ^ 



Jn{p) 



27r2 k^ + [C„(p)]2 ' 



(28) 



mpTUn 



-A 



I '^pnh 



m. 



p' + 



rridh 2 



pn 



rrid 



The wave functions ipdik) and (f)n{p) in the equations (26)-(28), which determine 
the model wave function of the triton, were found by solving analytically the 
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corresp ending two-body problems. Note that the p — n interaction is relatively 
weak resulting to rather small deuteron binding energy. This interaction can be well 
described using even a simple separable rank-1 potential. Another situation exists 
in the case of the n — d interaction that is stronger than the p — n interaction. In 
addition to the n — c? bound ground state (the triton), the n — d interaction supports 
also one virtual state (that is reflected in a rather small value of the n — d scattering 
length). In this connection the n — d interaction can be described by a separable 
potential of the rank no less than 2 or by a local potential. 

Here, in the case of the triton, both the wave functions ipdik) and (j)n{p) were 
determined by solving the corresponding two-body bound-state problems with the 
local Hulthen (H) potentiaP^, 

Ar) = -^;o[exp(gr)-l]-' , (29) 

employed to describe the proton-neutron interaction, f^(r), and the effective in- 
teraction between the neutron and the deuteron (as a structureless object), v^{p). 
(The radius vector variables r and p in configuration space correspond to the vari- 
ables k and p in momentum space.) 

4. Calculations and discussion of results 

The electric dipole polarizability of the triton q;£(^H) was calculated by the formulae 

(25) - (28). 

The parameters of the p — n interaction potentials were fitted to the experimental 
values of the deuteron binding energy and the triplet p — n scattering length ^ap„, 

= 2.224575(9) MeV (Ref. 27), ^a^^ = 5.424(3) fm (Ref. 28). (30) 

The parameters of the effective n — d interaction potential (29) were determined 
using the experimental values of the separation energy of the neutron (needed to 
remove one neutron from the triton, Bn = Bf- Bd) and the doublet n — d scattering 
length ^ttnd, 

Bt = 8.481855(13) MeV (Ref. 29) , ^Und - 0.65(4) fm (Ref. 30). (31) 

The fitted values of the parameters (7 = Voq~^ and /3 — q + k) oi the local 
Hulthen p — n and n — d ineraction potentials (29), Vpn and i'„,were found to be 

2/Xpn7pn = 1.3184 fm , /3pn = 1.3146 fm~^ ; 2//„d7„ = 5.4689 fm , ^„ = 0.9552 fm~^ . 

(32) 

For the deuteron, the electric dipole polarizability obtained from Eq.(3) for 
the local Hulthen potential with the parameters (32) takes the value q;^(^H) = 
0.6442 fm^. It is known '''■^^ that the S-wave asymptotic normalization constant 
AsC^lV) accounts for most of the polarizability q;e(^H) having regard to a rather 
high probability of that the slightly bound nucleons in the deuteron are at distances 
outside of the range of the nuclear force. If both the parameters of the potential are 
fitted to only the data that concerns to the bound p — n state, for example, with 
the use of the binding energy Bd (30) and 



AsCR) = 0.8845(8) fm'^/^ (Ref. 32) , 



(33) 
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we find for the deuteron polarizability the value 

a^CR) = 0.6292 fm^ . (34) 

With the use of the tensor separable potential that corresponds to ^4^ = 
0.8843 fm~^/^ the deuteron polarizability increases still further — to the value 
aE — 0.6311 fm^ (Ref. 8) — approaching the values obtained with the realistic 
potentials, = 0.6328(17) fm^ (Ref. 7). 

In this paper the electric dipole polarizability of the triton given by the formulae 
(25) ~ (28) has been calculated with the use of the local Hulthen potential that 
allows of finding the wave functions ipd{k) and 0„(p) in the analytical form. Fitting 
the parameters of the potential to the data for Bd and ^ap„ of the p — n system (30) 
and for Bn and "^and of the n — d system, (30) and (31), we have obtained for the 
triton polarizability the value 

a|(3H) = 0.235 fm^ . (35) 

In this case, the direct and exchange terms in (25) are found to contribute to the 
triton polarizability almost equally: 

a%{pnn) = 0.109 fm=^ , a^Hpnn) = 0.126 fm^ . 

If the parameters of the p — n interaction potential are fitted to the data 
characterizing only p — n bound state (the deuteron), Bd and AsC^H) (the S-wme 
asymptotic normalization) , we obtain a slightly less value of the triton polarizability, 
^^(^H) = 0.225 fm^ . The decrease of the polarizability in this case occurs due to 
the decrease of the quantity As (the experimental value of As (33) is less than the 
value of As for the Hulthen wave function corresponding to the potential fitted to 
the data Bd and ^Upn (^s(^H) = 0.8960 fm'^/^)). 

Unfortunately, the value ^^(^H) has not been measured yet. Although our ten- 
tative estimate of the triton polarizability, carried out on the basis of the direct 
calculation by the formula (21) but using the cluster wave function, leads to the 
result for ^^(^H) that is consistent with one of two data for the polarizability of the 
mirror nucleus ^He (Ref. 1), a new calculation of a£;('^H) with the use of a more 
reasonable triton wave function in (21) would be valuable. 



5. Summary and conclusions 

Leaning upon on the analytical structure of the three-body transition matrix, a 
consistent formalism for determination of the electric dipole polarizability of a three- 
hadron bound complex consisting of one charged and two neutral particles and 
having only one stable bound state has been worked out. A simple expression for 
the electric dipole polarizability of the three-hadron bound system has been derived 
assuming that the higher orbital components of the three-body off-shell transition 
matrix at the negative energy E — —Bq are negligibly small. In this case, the 
polarizability is expressed in terms of the first partial derivatives of the bound-state 
wave function with respect to the Jacobi momentum variables of the complex (Eq. 
(21)). 
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Applying the cluster model for the triton wave function and the Hulthen in- 
teraction potential to describe the p — n and n — d bound systems, we have 
obtained for the electric polarizability of the triton the approximate estimate 
q;^(^H) = 0.23 fm^(using the low energy data for Bd and ^ap„ of the p — n system 
and for Bn and '^and of the n — d system). Under conditions that there is presently 
no direct measurement of the quantity a£'(^H) and the results of the experiments for 
the polarizability of the mirror nucleus '^He (Refs. 3 and 4), are inconsistent, more 
complicated calculations of the triton polarizability on the basis of the expression 
(21) and modern wave functions for the bound state are worth to be performed. 
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